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Abstract 

This article presents an original code for Big Bang Nucleosynthesis in a baryon inhomogeneous 
model of the universe. In this code neutron diffusion between high and low baryon density regions 



{Sj . is calculated simultaneously with the nuclear reactions and weak decays that compose the nucle- 

>, 

osynthesis process. The size of the model determines the time when neutron diffusion becomes 
VO ' ' 

significant. This article describes in detail how the time of neutron diffusion relative to the time 

\q ; 

of nucleosynthesis affects the final abundances of 4 He, deuterium and 7 Li. These results will be 
compared with the most recent observational constraints of 4 He, deuterium and 7 Li. This inho- 
Qh' mogeneous model has He and deuterium constraints in concordance for baryon to photon ratio 

2 ■ n = (4.3 — 12.3) x 10 -10 7 Li constraints are brought into concordance with the other isotope con- 

straints by including a depletion factor as high as 5.9. These ranges for the baryon to photon ratio 
and for the depletion factor are larger than the ranges from a Standard Big Bang Nucleosynthesis 
$_i ■ model. 
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I. INTRODUCTION 



Big Bang Nucleosynthesis ( BBN ) is the primary mechanism of the creation of the 
lightest isotope species . At temperatures of the universe T m 100 GK baryonic matter 
consisted mostly of free neutrons and protons in thermal equilibrium with each other via 
weak interconversion reactions. Weak freeze-out occurs when the temperature falls to T px 
13 GK and the interconversion reactions fall out of equilibrium. Between T « 13 GK and 
T ph 0.9 GK only neutron decay changes the neutron and proton abundances. Then nuclear 
reactions become significant, forming heavier and heavier nuclei. Nearly all free neutrons 
at the time of nucleosynthesis are incorporated into 4 He nuclei because of the large binding 
energy of that nuclei. The amount of free neutrons at that time depends on the neutron 
lifetime r n . BBN is also the only source of deuterium production, and a significant source 
of 7 Li production. 

The nuclear reaction rates depend on the baryon energy density p&, equivalently the 
baryon to photon ratio rj. The abundance results of 4 He, deuterium and 7 Li can then be 
compared with measurements to put observational constraints on the value of rj. BBN 
constraints on rj can be compared with constraints derived from Cosmic Microwave Back- 



ground measurements 



Acoustic oscillations in the CMB angular power 
spectrum are fitted with spherical harmonic functions that depend on several cosmological 
parameters, including the density factor Q^h 2 . The most recent CMB measurements set 
Q B h 2 = 0.0224 ± 0.0009 3, corresponding to r] = (5.9 - 6.4) x 10^ 10 . 

The Standard Big Bang Nucleosynthesis ( SBBN ) model is the simplest BBN model. 
In SBBN all constituents are homogeneously and isotropically distributed. Parameters that 
define the SBBN model are rj, r n , and the number of neutrino species N v . But a variety of 
models alternative to SBBN can be fashioned by adding in other parameters. The ability for 
BBN to constrain the value of rj depends on the reliability of isotope measurements. Isotope 
observational constraints on rj have frequently appeared not to be in concordance with each 
other when applied to the SBBN model. The possibilit y of alternative BBN mode 



discrepencies has then been considered 
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24j,|25(. A good understanding of alternative models should then be maintained. 



Figure ^ shows graphs for the mass fraction X4 He of 4 He and the abundance ratios 
Y(d)/Y(p) and Y(jLi)/Y(p) of deuterium and 7 Li, all as functions of rj. These graphs 
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correspond to an SBBN model. The SBBN code used for Figure [T] has been used by this 
author in previous articles 

4 He is measured in low metallicity extragalactic HII regions. There is disagreement over 
how to extrapolate data points to zero metallicity. In some studies extrapolations have 
led to a higher mass fraction value of around 0.244 



value is a lower 0.234 
0.239±0.002 
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23, 123, 123, while in other studies the 
32I ] . The most recent measurements of X-i He have been a lower 
and a higher 0.242±0.002 Q. But the extent of systematic errors in these 
results is controversial. Olive et al 3^ have used a compromise value 0.238±0.005 combining 
both high and low measurements due the uncertainty in systematic error. Recently Olive 
and Skillman j^l try to quantify uncertainties due to systematic error, reporting a large 
range 0.232 < X4 Hc < 0.258. This range should eventually go down as the quantification of 
systematic errors improves. Figure ^ shows the 2a ranges by Luridinia et al j^] and Izotov 



and Thuan 



341 ] (IT04) combined, corresponding to a range of 77 = (2.2 — 6.1) x 10 



-10 



The deuterium measurement shown in Figure Q is the weighted mean of five Quasi-Stellar 
Objects ( QSO )'s done by Kirkman et al Q- This abundance^tio_T(d)/Y"^^ - n 7C+0 44 



10 5 is in good agreement with many previous measurements 



38, 39, 




40, 41 



— 2.78_o. 38 x 
42j . But Rugers 

and Hogan [43] measured Y(d)/Y(p) an order of magnitude greater, at (1.9 ± 0.4) x 10~ 4 . 
The abundance ratio by Kirkman et al corresponds to rj = (5.6 — 6.7) x 10~ 10 , which is in 
good agreement with the CMB results. 

Ryan et al [4^ measure 7 Li by looking at a group of very metal-poor stars and account- 
ing for various systematic errors to derive a value Y( 7 Li)/Y(p) = 1.23^032 x 10~ 10 . This 
measurement has a smaller magnitude and value of a than preceding measurements j3, 3 . 
The largest uncertainty in the calculation of this abundance range is the uncertainty in 



determining the effective temperature T e g of the stars. Melendez & Ramirez |47[ make new 
calculations of T e g- and get higher temperatures than Ryan et al for lower metallicity stars. 
Melendez & Ramirez then derive a larger value Y( 7 Li)/Y(p) = 2.34^Qg 4 x 10~ 10 , also with 
a larger 2a error. Figure Q shows the measurements by both Ryan et al and Melendez & 
Ramirez. Ryan et al's measurement corresponds to rj = (1.6 — 4.2) x 10~ 10 while Melen- 
dez & Ramirez's measurement can correspond to two ranges, rj = (1.1 — 2.0) x 10~ 10 and 
v = (3.3-6.0) x 10- 10 . 

This measurement of 4 He by IT04 is in concordance with the deuterium measurement of 
Kirkman et al only at its 2a range, for a narrow range rj = (5.6 — 6.1) x 10~ 10 . The 7 Li 
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constraints by Melendez & Ramirez is in concordance with the deuterium measurement also 
only at its 2a range, while the 7 Li constraints by Ryan et al have no region of concordance 
at all. A depletion factor from stellar evolution could improve concordance between the 
deuterium constraints and the 7 Li constraints by Melendez & Ramirez. A factor of 2.8 
would resolve the discrepency in the case of Ryan et al's constraints. But models for 7 Li 
depletion in stars and measurements of a depletion factor remain controversial. 

This article focuses on the particular alternative model of Big Bang Nucleosynthesis with 



an Inhomogeneous baryon distribution 
original code written by this author [4 



IBBN ). The IBBN code used in this article is an 
hereafter known as the Texas IBBN code. Upon 
publication of this article the code will be made publically available at the author's website 
49| . The Texas IBBN code can serve as a consistency check against other IBBN codes, and 
against SBBN codes as well when run in its small distance scale limit. 

Section |n] is a summary of the history of IBBN research, emphasizing developments that 
are significant to the way the Texas IBBN code is constructed. Section ITTT1 lists the specific 
details of the IBBN model used for this article. Section HVl shows the final abundance results 
of the IBBN code for a range of distance scale and baryon to photon ratio rj. This section 
discusses how the time of neutron diffusion relative to weak freeze-out and nucleosynthesis 
significantly affects the final isotope abundances the code produces. The description of this 
relation in this article is a useful guide for how baryonic matter flows and is processed in 
an IBBN model. In Section El the IBBN model will be compared with the most recent 
constraints on 4 He, deuterium and 7 Li. For certain IBBN parameter values the acceptable 
range of r\ from 4 He and deuterium constraints is widened. The IBBN model also permits a 
large range of Li depletion factor that is of particular interest. 

II. DEVELOPMENT OF THE IBBN CODE 

Various theories of baryogenesis lead to inhomogeneous distributions of free neutrons 
and protons by the time of nucleosynthesis. The distributions can be modelled with many 
different symmetries. Baryon inhomogeneities may arise from a first order quark hadron 
phase transition [sfll . Isil . Transport of baryon number between quark gluon phase and 
hadronic phase is inefficient, leading to concentration of baryon number in the last remaining 
regions regions of quark gluon plasma j&J. The magnitude of the bubble surface tension 
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determines if the quark gluon plasma regions form in centrallly condensed spherical bubbles 



13, 



or cylindrical filaments |55j. A cosmic string moving through matter during the 
quark hadron phase transition can also leave wakes of matter that remain in the quark 
gluon plasma phase longer than in the regions outside the wakes, forming sheets of planar 

571 ] . B aryon inhomogeneity may also form during the earlier 
A first order phase transition would proceed by 



symmetric inhomogeneity 
electroweak phase transition 
bubble nucleation. Particles in the plasma interact with the bubble walls in a CP violating 

□n 

matter, leading to a baryon asymmetry forming along the walls l60Ll61j ]. The baryon density 
is in the form of high density shells, spherical or cylindrical IfSfsl ] . Baryon inhomogeneities 
can also arise from phase transitions involving inflation-generated isocurvature fluctuations 



15711 . Baryoi 



62| . or kaon condensation phase 



64. 



651 . |66j treated regions of different 



The earliest articles on inhomogeneous codes 
density as separate SBBN models. They would run a model with a high value of pb, then 
a model with a low value, and then average the mass fractions from each model together, 
weighting each on how large each density region was. Applegate, Hogan and Scherrer 
considered the possibility of nucleons diffusing from high density regions to low density 
regions. Neutrons diffuse by scattering off of electrons and protons. Protons scatter off of 
neutrons and Coulomb scatter off of electrons, but the mean free path of protons is about 
10 6 times smaller than that for neutrons because of the Coulomb scattering. Diffusion of 
other isotopes is negligible compared to neutron scattering because the isotopes are much 
more massive. 

In early IBBN codes that featured neutron diffusion 0,0,0 the diffusion part is run 
first, at early times and high temperatures. Then nucleosynthesis within the regions is al- 
lowed to run. In their IBBN code Kurkio-Suonio et al (KMCRW88) made the significant 
innovation of having neutron diffusion occur both before and during nucleosynthesis. This 
code was for planar symmetric baryon inhomogeneity, and was split in a uniform grid of 
20 zones. Kurki-Suonio and Matzner 01 (KM89) and Kurki-Suonio et al j72] (KMOS90) 
looked at cylindrical and spherical models, using uniform grids as well. But for larger ratios 
between high and low densities, or lower volume fractions of high density region, the number 
of zones needed for the code to run accurately increased considerably. The codes used by 
Kurki-Suonio and Matzner (KM90) and Mathews et al 0, 75] instead use nonuniform 
grids, with a greater number of narrower zones around the boundary between high and low 



5 



density regions, where they are needed. Mathews et al 74] halve the width of a zone the 
closer the zone is to the boundary. KM90 use a stretching function to make a grid of 64 
zones that get very narrow around the boundary. 



III. THE IBBN MODEL 

In an IBBN model the universe is represented as a lattice of baryon inhomogeneous 
regions. An IBBN code models one region in that lattice. The inhomogeneity can have 
planar, cylindrical or spherical symmetry. The Texas IBBN code has been used to model 
condensed spheres jzf| and cylinders ( a high density core ) and spherical and cylindrical 



771 ] shells ( a high density outer layer ). The parameters that define an IBBN model are 
the baryon to photon ratio rj, the distance scale the density contrast R p , and the volume 
fraction f v . The distance scale is the initial size of the model at a chosen time. In this article 
that time is the starting time of a run, when the temperature T = 100 GK. The density 
contrast is the initial ratio of high baryon density to low density. The volume fraction is the 
fraction of the model occupied by the high density region. f v is parametrized such that it 
corresponds to a specific radius. 

A cylindrical shell model will be used in this article. This symmetry has been used by 
Orito et al 3] and Lara 2004 \r\ . The isotope abundance results are represented as contour 
maps in a parameter space defined by 7] and r^. The values of the remaining parameters are 



taken from Orito et al 



55]. 



R p = 10 6 
1 - \/\ - fv = 0.075 

The contour lines of abundance values to be discussed in Sections HYI and El are most greatly 
exaggerated in a cylindrical shell model with the parameter values from Orito et al, meaning 
that observational constraints will be satisfied for the highest possible values of fig/i 2 ( t] ). 
The parametrization of f v means that the thickness of the high density outer shell equals 
0.075 the radius of the whole model. For the neutron lifetime the most recent world average 
T n = 885.7 seconds 178] is used. 

The model is divided into a core and 63 cylindrical shells. These zones need to be thin at 
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the boundary radius r& between high and low density to accurately model neutron diffusion. 
The Texas IBBN code uses the stretching function from KM90 [73] to set the radii of the 
shells. 



= an) + 77- i 
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£(r) is the shell number out from the center, with a radius r in normalized units that range 
from to 64. The boundary radius r& = 59.2 as determined by the value of /„ Figure El 
shows how £ (r) maps onto r 

Appendix |X] describes in detail the method the calculations are made for each timestep 
in the run. 

IV. RESULTS 

Figures OH are contour maps of the overall mass fraction X4 He and abundance ratios 
Y(d)/Y(p) and Y( 7 Li)/Y(p) at the end of the Texas IBBN code's run, drawn in a parameter 
space defined by rj and r^. In Figure El the abundance ratios of both 7 Li and 7 Be are shown 
combined, as all the 7 Be has decayed to 7 Li by now. Neutron diffusion starts at the boundary 
between the high density outer region and the low density inner region, and then progresses 
outwards to the outermost shell and inwards to the core. The time neutron diffusion takes to 
homogenize neutrons determines the shapes of the contour lines shown in Figures The 
two milestone times in element synthesis are the times of weak freeze-out and nucleosynthesis. 
The contour lines can be described in terms of whether neutron diffusion occurs before weak 
freeze-out, between weak freeze-out and nucleosynthesis, or after nucleosynthesis. 

A. Before Weak Freeze- Out 

For the smallest distance scales neutron diffusion homogenizes neutrons very early in 
a run. Protons are still coupled with neutrons via the interconversion reactions. In the 
high density outers shells these interconversion reactions run in the direction of converting 
protons to neutrons, to keep up with neutron diffusion. The protons converted to neutrons 
diffuse to the low density inner shells, where the the interconversion reactions run in the 
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opposite direction, converting neutrons to protons. Protons are then homogenized along 
with neutrons, and the final abundances are the same as the abundances from an SBBN 
model. 

For larger neutron diffusion takes longer to affect all shells of the model. At a distance 
scale of around 1600 cm the time when diffusion ends coincides with weak freeze-out. Protons 
are not as coupled with neutrons as with smaller distance scales, and so are not completely 
homogenized by the time when neutrons have been homogenized. A larger proton density 
makes nucleosynthesis occur earlier in the outer shells. For a given value of r\ then the final 
abundance results are the results from an SBBN model with earlier nucleosynthesis: greater 
4 He, lesser deuterium, and greater 7 Li and 7 Be production. That corresponds to the shift in 
the contour lines of Figures EE to lower rj for distance scales from 1600 cm to 25000 cm. 

B. Between Weak Freeze-Out and Nucleosynthesis: 7 Li and 7 Be 

In models with r, from « 25000 cm to 3.2 xlO 5 cm neutrons are homogenized at a time 
in between weak freeze-out and nucleosynthesis. 

If rj ~ 25000 cm, neutron diffusion becomes significant everywhere right around the time 
of weak freeze-out. Figures EE show an example of how the various reactions in the code 
interact with one another. Figures EE correspond to shell number 62, a high density outer 
shell that is two shells away from the outer edge of the model. The reaction rates are 
normalized to the average baryon number density n&o and the expansion rate of the universe 
d R . 

In FigureEJthe neutron diffusion rates peak at around T = 10.0 GK and remain large up 
to T = 3.0 GK. The diffusion rate from shell 62 out to shell 61 is larger than the rate from 
shell 63 into shell 62 all through that time. The net effect is outflow of neutrons from shell 
62, as it is happening in all high density shells at this time. The peak temperature T = 
10.0 GK is just after the temperature T « 13 GK of weak freeze-out. Figure E] shows the 
rates for the reactions that convert neutrons to protons ( n — > p ) and the rates that convert 
protons to neutrons ( p — > n ) in short dashed lines. These rates are the same as they would 
be in the SBBN model. So no proton redistribution via these reactions is possible, and the 
proton number density in shell 62 remains high. 

Figure U\ shows the nuclear reaction rate n + p^d + 7in short dashed lines. This 
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reaction falls out of Nuclear Statistical Equilibrium ( NSE ) at T — 0.9 GK, starting off the 
chain of nucleosynthesis. Because the proton number density in the outer shells is high the 
nuclear reactions go at faster rates than they would in the SBBN model. Nucleosynthesis 
then occurs slightly earlier in the outer shells, depleting neutrons there. This deficit of 
neutrons leads to back diffusion. Figure [7| shows the rates of diffusion from shell 61 into 
shell 62, and from shell 62 out to shell 61. The net effect is now a concentration of neutrons 
in the high density shells. Nearly all nucleosynthesis is concentrated in the outer shells. 

7 Li is created primarily by the nuclear reactions t + 4 He <-> 7 Li + 7 and n + 7 Be «-> p + 
7 Li, and destroyed primarily by the reactions p + 7 Li <-> 2( 4 He ) and d + 7 Li <-> n + 2( 4 He 
). The depletion reaction p + 7 Li <-> 2( 4 He ) dominates over other reactions involving 7 Li. 
7 Be is created primarily by 3 He + 4 He <-> 7 Be + 7 and destroyed primarily by n + 7 Be <-> p + 
7 Li. In contrast to 7 Li the creation reaction of 7 Be dominates over the destruction reaction, 
and greater 4 He production in the high density shells magnifies the dominance even further. 
Figure |H1 shows the number densities of 7 Li and 7 Be as functions of radius. The number 
density of 7 Be is considerably larger in the high density outer shells than in the rest of the 
model. Due to this greater 7 Be production the contour lines in Figure El have a larger shift 
to lower rj than the contour lines in Figures EE1 

C. Between Weak Freeze-Out and Nucleosynthesis: 4 He and deuterium 

In models with rj from m 25000 cm to 10 5 cm the proton number density is unchanged 
from the time of weak freeze-out to nucleosynthesis, except for a slight increase due to 
neutron decay. For this range of rj the contour lines in Figures EK lie along nearly constant 
values of rj. 

In models with r« = 10 5 cm the amount of time needed for back diffusion to affect all 
shells is the same as the duration time of nucleosynthesis. For larger distance scales the shells 
furthest from the boundary are not as well coupled by back diffusion to the boundary shells. 
Nucleosynthesis becomes concentrated in the shells immediately around the boundary. This 
concentration leads to an overall drop in 4 He production. For from « 10 5 cm to 3.2 xlO 5 
cm the contour lines in Figures EE shift to higher rj. Figure El shows the final number density 
of 4 He as a function of radius for rj « 3.2 x 10 5 cm, with 4 He very concentrated around the 
boundary. 
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For models r« > 3.2 x 10 5 cm diffusion cannot homogenize neutrons before nucleosynthesis. 
A larger neutron number density remains in the outermost high density shells, and a lower 
density in the low density core and innermost shells. The larger neutron number density 
leads to greater 4 He production in the outermost shells. Figure E] shows the final number 
density of 4 He for = 2.0 x 10 6 cm. There is greater 4 He production around the boundary 
and the outermost shells and a trough of lower production in between. The overall 4 He mass 
fraction increases again, For r« > 3.2 x 10 5 cm the contour lines in Figure El and Figure El 
shift to lower T]. 

Decreased 4 He production tends to be accompanied by increased deuterium production. 
FigurelTUlshows the final number density of deuterium for r, pa 3.2 x 10 5 cm and = 2.0 x 10 6 
cm. In the radii corresponding to the trough of 4 He production in Figure El Figure El has 
a peak in deuterium production. The contour lines in Figure 0] shift to higher r\ for from 
pa 10 5 cm to 3.2 xlO 5 cm, just as in Figure El and Figure El But for from pa 3.2 x 10 5 cm 
to 2.0 xlO 6 cm the deuterium contour lines still shift to higher rj because of the increased 
deuterium production shown in Figure E3 

D. After Nucleosynthesis 

At rj pa 2.0 x 10 6 cm neutron diffusion peaks at the same time as nucleosynthesis. For 
models with larger r« neutron diffusion becomes less significant. More neutrons initially 
in the high density outer region remain there, increasing 4 He production. The trough in 
Figure El disappears and so deuterium production decreases. In Figure |U the deuterium 
contour lines shift to lower 77 to coincide with the contour shifts in Figure El and Figure El 
The largest models behave as two separate SBBN models; a high density SBBN model with 
considerable 4 He and 7 Li+ 7 Be production and minimal deuterium production, and a low 
density SBBN model with minimal 4 He and 7 Li+ 7 Be production and substantial deuterium 
production. Final results are the average results from the two models. 

E. Generalization 

The contour maps shown in Figures EH are for a specific IBBN model. If the model 
geometry is changed or if the values of the other parameters, the density contrast R p and 
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the volume fraction f v , are changed the shifts in the contour lines become more or less 
exaggerated. But the basic shapes of the contour lines persist. For all geometries and values 
of R p and /„ there will be a range of distance scale where neutron homogenization occurs 
in the interim between weak freeze-out and nucleosynthesis, leading to the shift to lower 77 
as shown in this article's model for r, pa 25000 cm. There will also be a range of where 
neutron diffusion coincides with nucleosynthesis. A trough of lower 4 He production between 
the boundary and the high density shells furthest from the boundary develops in this range, 
like the trough shown in Figure 01 IBBN models will then have a distance scale where the 
contour lines of 4 He and deuterium diverse. For a talk at the Sixth ResCEU International 
Symposium |77( this author looked at models with the geometries of condensed cylinders, 
condensed spheres, and spherical shel, as we.1 as cylindrical she.., The va.nes of R e and /, 
used by Orito et al 55J were used in those runs. The contour maps in all the models showed 
the same features as seen in Figures EE 

V. OBSERVATIONAL CONSTRAINTS 

Figure HP shows the observational constraints from Figure U a ppli ed to the contour maps 
of Figures OK The maximum X4 He < 0.246 constraint from IT04 |34| and the 7 Li constraints 
from Ryan et al [44] are shown in Figure ITT1 

Regions of concordance between the IT04 4 He maximum constraint and the deuterium 
constraints are shown in yellow. A concordance region exists for distance scales < 5000 
cm and rj = (5.6 — 6.1) x 10~ 10 . These limits on 77 are the same limits as seen in SBBN 
models. The maximum limit of r« is set by the shift to lower 77 as neutron diffusion occurs 
closer to weak freeze-out. The X4 He = 0.246 contour have a greater shift than the contour 
lines for the deuterium constraints, because of increased 4 He production in the outer region. 

Another region of concordance appears for rj = (1.3 — 6.0) x 10 5 cm, when the contour 
lines shift to higher 77 due to the concentration of nucleosynthesis along the boundary. The 
upper cutoff of is determined by the condition when a trough as shown in Figure El exists 
in the 4 He abundance distribution. Greater 4 He production in the outermost shells cause the 
A4 He = 0.246 contour to shift to lower 77 while greater deuterium production in the trough 
cause the deuterium contour lines to remain shifted to higher 7/. The acceptable range of 77 
is (4.3 - 12.0) x 10~ 10 , larger than in the SBBN case. 
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The 7 Li constraints from Ryan et al Q] are shown in darkest green in Figure ^2 The 
contour lines for 7 Li tend to shift in the same direction with the contour lines of 4 He and 
deuterium. So the 7 Li constraints do not have a region of concordance with the 4 He and 
deuterium constraints for this IBBN model, and the lack of a region of concordance persists 
for other geometries and parameter values. Figure ^2 also shows the region of the 7 Li con- 
straints with a depletion factor of 2.8. That depletion factor would bring the 7 Li constraints 
in concordance with the other isotopes for distances scales T{ < 5000 cm. For the region of 
concordance corresponding to 7*j = (1.3 — 6.0) x 10 5 cm a larger depletion factor of 5.9 is 
needed. The greater production of 7 Be shown in Figure |H1 leads to the larger shift to lower 
7] in the 7 Li contour lines compared to the 4 He and deuterium contour lines, and the larger 
depletion factor. Figure ^2 shows the region of the 7 Li constraints with the depletion factor 
of 5.9. The benefit of IBBN models then is to allow for a larger range of 7 Li depletion factor 
than permitted by the SBBN model. 

Figure ITT1 is similar to Figure 2 from the proceedings article Lara 2004 \r\ . Differences 



between the figures include use of the newer Xi He < 0.246 constraint [34] in place of the 
Xi-Re < 0.248 constraint [35! . The method of calculating the diffusion coefficients j^J is also 
newer than the method [3] used in Lara 2004. The neutron lifetime r n = 885.7 seconds [7^ 
was also updated for this article. 

Figure IT2l shows the 7 Li constraints from Melendez & Ramirez 4JJ. The outermost edge 
of these la constraints has concordance with most of the concordance region between 4 He 
and deuterium for < 5000 cm. With a small depletion factor of 1.35 these 7 Li constraints 
cover the whole region of concordance. A larger depletion factor of 2.8 is needed to cover 
the region of concordance corresponding to 7*j = (1.3 — 6.0) x 10 5 cm. The regions of concor- 
dance between 4 He and deuteurium are controversial because of considerable disagreement 
regarding 4 He constraints. Nonetheless both Figure E] and Figure ^] show that Inhomoge- 
neous Big Bang Nucleosynthesis allows for a larger range of acceptable 7 Li depletion factor 
to bring deuterium and 7 Li in concordance with each other, due to the greater shift in 7 Li 
contour lines to lower 77 for distance scales Ti from m 1600 cm to 10 5 cm. 
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VI. CONCLUSIONS 



The Texas IBBN code is an original code written such that the weak and nuclear reactions 
of element synthesis are coupled with neutron diffusion. The time of neutron diffusion 
relative to the times of weak freeze-out and nucleosynthesis have a significant influence on 
the final production amounts of 4 He, deuterium, and 7 Li. Because diffusion is coupled to the 
reaction network the code correctly accounts for neutron back diffusion, wherein neutrons 
flow back into regions with higher proton density due to earlier nucleosynthesis in those 
regions. Back diffusion has an influence over the results especially when the time of neutron 
diffusion is close to the time of nucleosynthesis. Of most interest in the results is the larger 
range of depletion factor for 7 Li that the IBBN model permits over the SBBN model. 

In models where diffusion homogenizes the neutron distribution before weak freeze-out 
protons are coupled with the neutrons via the weak interconversion reactions. Protons are 
then redistributed. Proton redistribution is less effective in models with the time of diffusion 
closer to the time of weak freeze-out, leaving a higher proton density in the outer shells. 
Increasing proton density leads to earlier nucleosynthesis in the outer shells. Neutrons then 
back diffuse into the outer shells, concentrating nucleosynthesis there. Nucleosynthesis in 
the high density shells produces decreasing amounts of deuterium and increasing amounts of 
4 He and especially 7 Be. The increased production of 7 Be is significant in the determination 
of the depletion factor of 7 Li. 

For models with the time of diffusion close to the time of nucleosynthesis neutron back 
diffusion becomes less effective. Nucleosynthesis is concentrated in the volume fimmediately 
around the boundary. This concentration leads to decreasing 4 He and 7 Li+ 7 Be produc- 
tion, and increasing deuterium production. In models where the time of neutron diffusion 
coincides with nucleosynthesis neutrons are not homogenized during nucleosynthesis. An 
increasing neutron number density remains in the outermost shells as well as a decreasing 
number density in the innermost shells. 4 He, 7 Li and 7 Be production jumps in the high 
density outermost shells, and overall production of these isotopes increases. But between 
the boundary and the outermost shells are shells with a trough of low 4 He production. Deu- 
terium is produced in large amounts in that trough. The deuterium contour lines in Figure 0] 
diverge from the contour lines in Figures El and EJ 

For models where neutron diffusion peaks at the same time as nucleosynthesis the trough 
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in 4 He production has disappeared. Deuterium production decreases and the deuterium 
contour lines in Figure 0] are in line with the lines in Figures HU and EJ The divergence in 
the directions of contour lines is significant in setting constraints on 77 and in the IBBN 
model. 



Application of observational constraints to this 



BBN model found slivers of concordance 



and He constraints by IT04 



3- 



Con- 



between the most recent deuterium constraints 
cordance occurs for rj = (5.6 — 6.1) x 10~ 10 and fj < 5000 cm, and for rj = (4.3 — 12.3) x 10~ 10 
and r, = (1.3 — 6.0) x 10 5 cm. The point of divergence between the 4 He and deuterium con- 
tour lines sets the maximum limit of acceptable 77. The reliability of 4 He constraints remains 



controversial 



36|. 



Contour lines between 4 He, deuterium, and 7 Li run roughly par allel to each other. The 
region Figure E] marked by the 7 Li constraints by Ryan et al 44j then does not have an 
overlap with the slivers of concordance of 4 He and deuterium. A depletion factor of 2.8 
would bring concordance in both the cases of SBBN and the first region of concordance. 
But because of the larger shift of the 7 Li contour lines to lower 7] a larger depletion factor of 
5.9 is needed to bring the 7 Li constraints in agreement with the second region of concordance. 
Recent 7 Li constraints by Melendez & Ramirez (^tJ have weak concordance with 4 He and 
deuterium constraints in the SBBN case. But an IBBN model still allows for a larger range 
of depletion factor, up to 2.8, to have 7 Li be in concordance with 4 He and deuterium. 
The IBBN abundance results for 7 Li will be compared with new measurements of the 7 Li 

jrimordial abundance derived from the ratio ( 7 Li/ 6 Li) measured in the InterStellar Medium 
A new neutron lifetime r n = 878.5 ± 0.7± 0.3 seconds has recently been measured 

82 1 . Constraints on rj in an SBBN model have been reassessed with the new lifetime jssj ] . 
and the constraints on r] and r ?i in Figures ^2 and ^] will also be reassessed with the new 
lifetime in an upcoming article [84j. Additionally, this article will be followed up by articles 



4jj to both SBBN and IBBN 



applying an original solution of the neutrino heating effect 
models. 

VII. ACKNOWLEDGEMENTS 

This work was partially funded by National Science Foundation grants PHY 9800725, 
PHY 0102204, and PHY 035482. This author thanks the Center for Relativity at the 



14 



University of Texas at Austin for the opportunity to work on this research, and Professor 
Richard Matzner in particular for his help in preparing this article. This author also thanks 
Professor Toshitaka Kajino of the National Astronomical Observatory of Japan for his advice 
on what to focus in this article. 



APPENDIX A: TRACE OF THE TEXAS IBBN CODE 

The stretching function that sets the radii of the zones in the cylindrical shell model 



Ci 



1 



1 




arctan 



(r - r b ) 




+ 



r — r& 
C1C2 



(Al) 



was used by KM90 [73]. The radius r is normalized to equal 64 for the full radius of the 
model. Eq. IA1I maps radii r of the zone boundaries to unit values of £. Zones near have 
a width around C\. Zones far from r& have widths determined by C3 and a rate of zone-size 
change controlled by CV r& always corresponds to a unit value of £. For the model in 
this article there are 20 zones covering the high density outer shell and 44 covering the low 
density inner region. The baryon number densities n^-high i n the outer shell's zones and 
rib-iow in the inner region's zones are set 



nb ~ low ~ f v R P + (l-f v ) { } 

n b-high = Rp n b-low (A3) 

such that the number density averages out to n&o over the whole model. 
At any given timestep the code solves the differential equation [3] 



for the number density n(i, s) of isotope % in zone s. The first two terms correspond to the 
weak and nuclear reactions that destroy ( [ij] ) or create ( [kl] ) isotope i within zone s. 
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rib(s) is the total baryon number density in zone s and Y(i, s) is the abundance Y(i, s) of 
isotope i in zone s. 



The 3oiRn(i, s) term corresponds for the expansion of the universe, where R is the expansion 
coefficient of the universe and an = \nR. This term can be eliminated by transforming 
to comoving coordinates. From here on r will be in comoving coordinates. The last term 
corresponds to diffusion of isotope % into and out of zone s. The factor p depends on the 
geometry of the model, p = for planar symmetry, 1 for cylindrical symmetry and 2 for 
spherical symmetry. Currently only neutrons can diffuse in the Texas IBBN code. The 
neutron diffusion coefficient D n is calculated from the coefficients D ne for neutron electron 
scattering and D np for neutron proton scattering. 



1 



(A5) 



Banerjee and Chitre derived a master equation for the diffusion coefficient between two 
particles scattering off of each other, based on the first order Chapman-Enskog approxima- 
tion S6|. Kurki-Suonio et al (KAGMBCS92) [73], and Jedamzik and Rehm 24 1 derive the 
same equation for the diffusion coefficient D ne for neutron-electron scattering. 



3 pK c K 2 (z) n n 

D ne = ~\ ~ J=r= pr(l A6 

8 V 2 n e a nc y/zK 5/2 {z) n t 

K 2 (z) and K 5 / 2 (z) are modified Bessel functions of order 2 and 5/2, a ne is the transport 
cross section of the scattering and z = m e /kT. n n /n t is the neutron fraction of the total 
number of ALL particles. This fraction is of the order 10~ 10 and so can be ignored. For 
neutron-proton scattering Jedamzik and Rehm derive an updated expression for the 
diffusion coefficient D np 



16 



D 



3 c 1 / k B T 1 



(A7) 



I{a, b) 




up 



rriN is the nucleon mass. The parameters a s = -23.71 fm, r s = 2.73 fm, at = 5.432 fm, and 
r t = 1.749 fm come from singlet and triplet scattering. 

The Texas IBBN code progresses a timestep At m for each step m. Eq. IA4I is evolved 
using a implicit second order Runge-Kutta method |87J • To use this method, Eq. IA4I has to 
be linearized. The weak-nuclear reaction terms can be linearized in a manner similar to the 
linearization of abundances Y used by Wagoner 0] and this author 



At m _! is the time difference between step m — 1 and step m. For the diffusion term the zones 
are defined on a grid whose points r(s) correspond to the outer radii of zones s. Number 
densities n(i, s) are considered the number densities at the midpoint radius between the inner 
and outer radii of zone s. The points r(s) correspond to points in £(s) space a distance of 
one unit between each other. The first space derivative in the diffusion term in Eq. IA4I can 
be discretized as: 





[N i Y^- 1 (i,s)Y^(j,s)n mA (i,s) + 



(k, s)Y^ 1 (I, s)n mA (k, s) + 



(A8) 
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Note that the coefficient [r p D(dl;)/(dr)] depends on r. The (l/r p )(d/dr) can be rewritten 
as a partial derivative of r p+1 . One can then write the discretization of the second space 
derivative as: 



dn 
~dt 



d 



d(rP +1 ) 



r p D — 

or 



n[r(s)] — n[r(s — 1)] 



8-1/2 



dn 



~n7 = •■■(p + 1) 



(r*I>g),{n[r(« + |)]-n[r( S -|)]} 



(^gj^W^-^-n^ -§)] }' 
rP +1 (s) - rP +1 (s - 1) 

For an isotope i the densities ra[r(s + 1/2)], n[r(s — 1/2)] and n[r(s — 3/2)] are defined as 
n(i, s + 1), s) and n(z, s — 1) respectively. The coefficient [r 2 D(d^)/(dr)] s is calculated at 
r(s). One can apply this discretization to an implicit version of the diffision part of Eq. IA41 



n mA {hs) -ra m _i(i,s) (p+1) (r p £>f ) s 

At m _! ' ' ' r p +\s) - r p +\s - i) HmA[h S + L) 

(r p Df) + {r p Df) 
+(P + *) r^{s) S -rv+\s-T) n - A(i ' S) 

-n mA {i,s-l) (A9) 



rP +1 (s) - r p+l (s 



Any baryons that flow out beyond the distance scale are assumed to be replenished by 
baryons flowing in from other sets of shells, and r(0) = is the center of the shells. The 
code uses reflective boundary conditions at the endpoints of the grid. 
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n m A(i,l)-n m -i(i,l) (p + 1) {r p D§ ) 



-n mA (i, 1) 



(P + 1 )(^ 0P ;i 



rP +i(l) _ r p+i( ) 

^^2^64) _ (p+l)(r^§) 64 

At m _ x ~ " ' + rP+!(64) - rP^(63) nmA[h > 

(P+l) (ggf)g „ f . 6 ^ 
r^(64) - ^(63) m M ' J 

where r(64) = distance scale rj. The above equations can be applied to the diffusion of any 

isotope i, but only neutrons ( % = 1 ) diffuse for the results of this article. 

Eq. IA8I and Eq. IA9I combined together can be rewritten as a matrix equation for a new 

number density value n m A(i, s). The matrix consists of a 68 x 68 matrix for each of the 64 

zones, built from the terms in Eq. IA8I From Eq. IA9I come terms that couple n(l, s) with 

n(l, s + and n(l, s — 1) due to neutron diffusion. n mA (i, s) is then used in the following 

equation 



n m (i, s) = n m (i, s) + 



n mA {i, s) - n m -i(i, s) 



At r 



At m -i 

to calculate an interim value n m (i,s) of the number densities. This is the first step of 
the Runge-Kutta method, with h m (i,s) the first estimate of the values of n m (i,s) at time 
t m + At m . Using Y m (i, s) = n m (i, s)/n b (s) the code solves a second matrix equation 



At m 



n mB {i, S )-n m {i,s) _ y N M [yy^-lf. ,)yN j(q s r. , , 

Jv i y^(i |S )y T ^- 1 o-,s)n mfl O'.«)] 

■[iV fe F^- 1 (fc, S )F^(/, S )n mB (fc, S ) + 



N iYm k (k, s)^ 1-1 ^, s)n mB (l, s)] - 
(p + 1) (r p Df ) 

+(P + !) rP+ i (g) _ rP+ i (g _ i) s ) 

(p+i)^!)^ 

-n mB (i, s-1) 



rP +1 (s) - rP +1 (s - 1) 
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for new number density values n m B(i,s). At m is the time difference between step m and 
step m + 1. Final new values for n m+ i(«, s) at timestep m + 1 can then be calculated. 



n m +i(i, s) = n m (i, s) + - 



n m A(i, s) - n m _i(i,s) n mg (i, s) - n m (i, s) 



At™ (A10) 



At m _ 1 At m 

This is the second step ( "5" ) of the Runge-Kutta method. At the same time as with 
n m (i, s) the Texas IBBN evolves lni? and the electromagnetic plasma energy density p c+7 



d(kiR) 18 ~ . 
dt = A/ ^G(p, + p c + Pu ) 

dp e+1 



dt 



R R . 
R Pj ~ 3 r^ Pc + pc 



(All) 
(A12) 



also by the Runge-Kutta method. 

After both Runge-Kutta steps have been done the code determines the new baryon num- 
ber density n&(s) of each zone using 



68 



n b 



where Aj, is the atomic weight of isotope i. From rib(s) and n(i, s) the code can calculate 
Y(i,s). At any given time the abundance Y av (i) and mass fraction Xj of isotope i in the 
entire model can be calculated from Y(i, s) using 



Y av (i) 



j:T=iMs)[r p+1 (s)-rP+\s-l)] 
Aj = AjY av {i) 

These overall abundances and mass fractions can be shown as contour maps of the IBBN 
code's parameters, and compared to observational constraints. 
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FIG. 1: X 4Hc , Y(d)/YM and Y{ 7 U)/Y( V ) for the SBBN case. The graph for Xa Uc includes the 
measurements by IT04 ( 34] solid lines ) and by Luridinia et al ( 33] dashed lines ). The graph for 



Y 



7 Li)/y(p) shows measurements by both Ryan et al ( 



( 



441 ] solid lines ) and Melendez & Ramirez 



471 ] dashed lines ). 
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FIG. 2: The stretching function will map integer values of £ to values of radius r such that the 
r's near the boundary radius r&. are closely spaced. Here r& = 59.2. 
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V 

FIG. 3: The mass fraction Xa^ c in the IBBN code. The horizontal axis is for baryon-to-photon 
ratio i] and the vertical axis is for distance scale r\ in centimeters at temperature T = 100 GK. 
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FIG. 4: The log of abundance Y(d)/Y(p) 



28 




29 




FIG. 6: Neutron diffusion vs. the interconversion reactions for rj « 25000 cm. The rate of diffusion 
from shell 63 to shell 62 is shown in solid while the greater rate from shell 62 to shell 61 is shown 
in long dashed lines. The net result of diffusion is neutron depletion in shell 62. The forward and 
reverse direction rates of the neutron to proton conversion reactions are shown in short dashed 
lines. Neutrons and protons are decoupled during the peak time of diffusion. 
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FIG. 7: Back diffusion for r\ m 25000 cm. The rate of diffusion from shell 62 to shell 63 is shown 
in solid while the greater rate from shell 61 to shell 62 is shown in long dashed lines. The direction 
of diffusion is now reversed, resulting in a net increase of neutrons in shell 62. Back diffusion 
coincides with the time when the nuclear reaction n + p <-> d + 7 ( shown in short dashed lines ) 
falls out of NSE. 
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FIG. 8: The final number densities of Li and 7 Be produced in each shell of the model, for vi ~ 
25000 cm. 
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FIG. 9: The final number density of 4 He produced in each shell of the model. For r{ ~ 3.2 x 10 5 
cm 4 He production is concentrated surrounding the boundary. For r j = 2 x 10 6 cm 4 He production 
has incresed in the outermost shells due to neutrons not homogenizing. 
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FIG. 10: The final number density of deuterium produced in each shell of the model, for n ~ 
3.2 x 10 5 cm and r j = 2 x 10 6 cm. Deuterium production remains considerably large in the trough 
of 4 He production that arises in Figure El 
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FIG. 11: ( Color Online ) Observational constraints of 4 He, deuterium, and 7 Li are shown on 
the IBBN cylindrical shell model with 1 — y/1 — f v = 0.075, and R p = 10 6 . The constraints 
y( 7 Li)/y(p) = 1.23^032 x 10~ 10 44] are shown. Also shown are 7 Li constraints with depletion 
factors of 2.8 and 5.9. 
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FIG. 12: ( Color Online ) Same as in Figure HT| but showing the constraints Y( 7 Li)/y(p) = 
2.341 

a 96 x I*-* 10 m3 • Also shown are these constraints with depletion factors 1.35 and 2.8. 



36 



